Simplified approach to double jumps for fluorescing dipole-dipole interacting atoms 
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A simplified sciieme for tiie investigation of cooperative effects in the quantum jump statistics of small num- 
bers of fluorescing atoms and ions in a trap is presented. It allows the analytic treatment of three dipole-dipole 
interacting four-level systems which model the relevant level scheme of Ba"*' ions. For the latter, a huge rate 
of double and triple jumps was reported in a former experiment and the huge rate was attributed to the dipole- 
dipole interaction. Our theoretical results show that the effect of the dipole-dipole interaction on these rates is 
at most 5% and that for the parameter values of the experiment there is practically no effect. Consequently it 
seems that the dipole-dipole interaction can be ruled out as a possible explanation for the huge rates reported in 
the experiment. 

PACS numbers: 42.50.Ct, 42.50.Ar, 42.50.Fx 



I. INTRODUCTION 

The dipole-dipole interaction between atoms and molecules 
is of fundamental importance in nature as it gives rise to the 
all pervading van der Waals force. In physics, cooperative ef- 
fects in the radiative behaviour of atoms due to their mutual 
dipole-dipole interaction have also attracted considerable in- 
terest in the literature 1 1 ], and they may play a role for possible 
quantum computers based on trapped ions or atoms. Atoms 
exhibiting macroscopic light and dark periods in their fluo- 
rescence may provide a sensitive test for such cooperative ef- 
fects. Such macroscopic light and dark periods can occur in 
a multi-level system if the electron is essentially shelved in a 
metastable state, thereby causing the photon emission to cease 
01. Two such systems accordingly exhibit a dark period, a 
bright period of the same intensity as that of a single system, 
and a bright period of double intensity. Three systems exhibit 
an additional bright period of threefold intensity. The dipole- 
dipole interaction may now alter the statistics of these periods. 

In an experiment with two and three Ba+ ions |3, 4] a 
large number of double and triple jumps, i.e. jumps by two 
or three intensity steps within a short resolution time, had 
been observed, by far exceeding the number expected for in- 
dependent atoms. Theoretically, the quantitative explanation 
of such large cooperative effects for distances of the order 
of tenjvave lengths of the strong transition proved difficult 
|5, 6, 7, 8, 9, 10]. On the other hand, experiments with dif- 
ferent ions showed no observable cooperative effects 1 1 1], in 
particular none were seen for Hg"*" for a distance of about 15 
wave lengths fl^. More recently, effects similar to reference 
[3J were found in an experiment with Ca+ ions 1 13], in con- 
trast to another, comparable, experiment jQ]. Neither were 
cooperative effects found experimentally in an extensive anal- 
ysis of the quantum jump statistics of two trapped Sr+ ions 
iTal . Skornia et al. 1 16] recently put forward a new proposal 
for observing the dipole-dipole interaction of two V systems. 

For two V systems numerical (ItI] and analytical |JJ1] inves- 
tigations of the effect of the dipole-dipole interaction showed 
an increase of up to 30 % in the double jump rate when com- 
pared to independent systems. However, the s^tems used in 
the experimental setups of references |3, Il3 11911 were not V 
systems so that a direct comparison between theory and ex- 



periment was not possible. For this reason the present au- 
thors extended their investigation to two other systems L2C|, 
namely a D shaped system modeling the Hg+ ions used in 
reference and a four-level system (see Figure ^ model- 
ing the Ba+ ions of references |3, 4]. For two D systems co- 
operative effects in the same order of magnitude as for the V 
systems were found for ion-distances of a few wavelengths of 
the laser-driven transition. For larger distances practically no 
effects where found, in agreement with the experiments 1 1^ 
and with the results of reference |21]. In contrast, only negli- 
gible effects for a wide range of ion-distances were found for 
two of the four level-systems. Although this result contradicts 
the findings of references ^ a dkect quantitative compar- 
ison with the experiments was not possible since explicit ex- 
perimental data were only provided for three Ba+ ions. For 
this reason three of the D and V systems were investigated 
in reference |22]. In comparison to two of either systems the 
cooperative effects found in this case are considerably higher, 
namely up to 170% deviation from the case of independent 
atoms. However, since the complexity increases dramatically 
for higher-level systems, this approach could not be applied to 
three of the four-level systems which we use to describe the 
situation of reference j3|] . 

In the present paper a simplified approach for the calcula- 
tion of the transition rate will be presented with which three 
four-level systems can now be treated analytically. This ap- 
proach is valid for atoms with a level structure in which the 
transitions between the different intensity periods take place 
incoherently, i. e. via decay or via incoherent driving. The 
transition rates for three dipole-interacting four-level systems 
will be calculated. Cooperative effects for this system are 
found to be less than 5% and negligible for the experimental 
parameters of reference (|3|. Consequently it seems that the 
dipole-dipole interaction can be ruled out as a possible expla- 
nation for the huge effects measured in the latter experiment. 

In section n] the Bloch equation approach is recapitulated. 
On this basis the new method is presented in section Hill and 
applied to the four-level systems in II VI In section|V]the pos- 
sibility of a translation of this method to V system type level 
structures is discussed. 
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FIG. 1: (a) Relevant level scheme of Ba+ For the effective 

four level system the circled levels are merged to a single level, (b) 
Effective four-level system for Ba+. Strong coherent driving of the 
1 1) — 13) transition by a laser, weak incoherent driving of the 1 1) — |4) 
transition by a lamp, weak decay of level 1 2) . 



n. BLOCK EQUATION APPROACH 

The fluorescence, i.e. the stochastic sequence of photon 
emissions, of a system consisting of a number of atoms with 
macroscopic bright and dark periods can be described by a 
telegraph process. This process is characterized by the tran- 
sition rates between the different intensity periods. In refer- 
ences fTsl l22il they were calculated for different model 
level systems and different numbers of atoms using a perturba- 
tion approach based on the Bloch equation of the correspond- 
ing systems. This approach will be illustrated in the following 
by applying it to the simple case of a single three level system 
in a D-type configuration as depicted in figure|2] 

The Bloch equations can be written in the compact form 

il 



HcondP - 



cond 



(1) 



where i/cond is the conditional Hamiltonian of the quantum 
jump approach |24], for this system given by 



and <^(p) is the the reset state, 

^(p) = Ai|l>(2|H2)(l|+A2|2)(3|p|3)(2|+^3|l)(3|H3)(l|. 

(3) 

The Rabi frequency fl^ and the Einstein coefficients 
^1,^2, ^3 are subject to the condition 



n3,A3 > ^1,^2 



(4) 



A detuning of the laser has been neglected for simplicity. If 
the small optical parameters Ai, A2 are neglected the system 
splits into independent subspaces. They are given by 



^o = {|2)}, 



^i = {|l>,|3)} 



(5) 



These subspaces can be associated with the periods of in- 
tensity in the sense that in a period li the system is mostly 
in the subspace S^i. Taking a state po,i in one of the subspaces 
S^i at a time to we calculate the state at a time tg + At later 
in perturbation theory with respect to the small parameters. 
The time interval At used here should be long compared to 
the mean time between the emission of two photons but short 
compared to the length of the intensity periods, 



A^^,n^^ < Ai < A^\A^^ 



(6) 



For the calculation the Bloch equation is written in a Liouvil- 
lean form. 



p^^p^{^aiA3,n3)+J^i{AuA2)}p. 
The density matrix at time to + At is then given by d 

p{to + At, po,j) = Pss.j + / dr e^«''^ipss,» , 



(7) 



(8) 



where pss.i is the quasi-steady state in subsystem S^i, i.e. a 
steady state of One can write 



1 



P, 



(9) 



with p containing the contributions from the eigenstates of ^0 
for non-zero eigenvalues. This leads to i20ll 

1 

p{t+At,pi)^Ps,,,+j2^^3Ps^.j^t+i<^~-^oy^p- (10) 



1 3) 



strong laser, ^3 



A. 



A. 



i?co„d = |[(A2+A3)|3)(3|+Ai|2)(2|]+^[|l)(3| + |3)(l|] 

(2) 



FIG. 2: Three-level system in D configuration with fast transitions 
(solid lines) and slow transitions (dashed lines). 
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The last term can be neglected and the coefficient 0;^ can 
therefore be interpreted as transition rate pij from intensity 
period li to period Ij. They can be calculated by means of the 
dual eigenstates for eigenvalue of ^0 |20|. For a single D 
system the quasi-steady states are given by 



Pss,o-|2)(2|, 
1 



[iAl + nl)\i){i\ + nl\3)(3\ 



i^3f^3|l)(3|-iA3f73|3)(l|] 



(11) 



(12) 



for the dark and the light period respectively. The correspond- 
ing dual states are 

p^,-|2)(2|, and = |1)(1| + |3)(3| (13) 

From il 1> and (I12> one finds 

^iPss,o = -Ai|2)(2|+Ai|l)(l| (14a) 

and 



iAz ^3^3 
2 ~t~ 2i^2 



y4.2 ~f" 2f^2 



(|1>(3|-|3)(1|). 



The transition rates are then calculated from 

Pi] = aij = Tr(p^j'Jfipss,i)- 

as 

Poi = ctQi = Ai 

and 



PW = "10 



A2nl 



(14b) 



(15) 



(16) 



(17) 



in agreement with the direct calculation of the transition rates 
via the quantum jump approach. 



III. NEW SIMPLIFIED APPROACH 

Due to the increased number of levels involved, a calcu- 
lation of the transition rates for three dipole-interacting four- 
level systems would, although in principal feasible with the 
methods introduced above, be even more laborious than for 
three three-level systems. It is, however, possible to read off 
the transition rates without having to carry out the full calcu- 
lation. One only needs the quasi-steady states of the corre- 
sponding subsystems. In following this simpler approach will 
be presented. 

By looking at equation (I14> one realizes that the last step 
in the calculation, namely the projection onto the dual eigen- 
states, although formally more satisfactory, was actually not 
necessary in order to gain the final result. The transition rates 



are already present as prefactors for some of the density ma- 
trix elements. In fact, ^1 can be interpreted as a transition 
operator Applying it to some state of the system yields the 
density matrix elements which are modified by the weak de- 
cays multiplied by the corresponding decay rates. They are 
positive for density matrix element which gain population and 
negative for those which loose population due to the decay. In 
the case in which one started with pss.o = |2) (2| one therefore 
has a term — ^i|2)(2|, which accounts for the loss of popu- 
lation of level 1 2), and a term Ai 1 1) (1 1 for the corresponding 
gain of population in the ground state. When starting with 
Pss,i the Einstein coefficient A2 for the decay from |3) to |2) 

has an additional factor ^^^^^^ for the quasi-steady state pop- 
ulation of level 1 3). The last two terms in equation il4h\ are 
due to the decay of the coherences between |1) and |3). 

From these considerations one is lead to a simple scheme 
for the evaluation of the transition rates. First one has to iden- 
tify the different independent subspaces for vanishing weak 
decay rates and calculate the quasi-steady states in these sub- 
spaces as in the above Bloch equation approach. For a single 
D system these are the states pss.o and pss,! for the subsystems 
associated with the dark and bright period, respectively. By 
looking at the level scheme one can then determine the pos- 
sible decay channels between the subsystems. In the present 
case this is a decay by A2 from |3) to |2) and a decay by Ai 
from 1 2) to The transition rates are then given by these 
decay rates multiplied with the steady state population of the 
decaying level. 

Physically this is quite intuitive: The transition rates are 
given by the corresponding decay rates multiplied by the mean 
occupation probabilities of the levels involved. 

The question is now if this approach can be extended to 
more complicated systems, especially to dipole-interacting D 
systems and to the four-level system for the description of 
Ba"*". This is indeed possible. For two dipole-interacting D 
systems for example the possible decays can be read off Fig- 
ure |3l which shows the level scheme in the Dicke basis given 
by 



|5) = |1)|1), |e2) = |2)|2), |e3) 
, , 1 



|3)|3) 



V2 
1 



(18) 



The easiest case is the transition rate poi for a transition from 
a dark period to a period of intensity /i. Here the relevant 
transitions are from |e2) to |si2) and |ai2). The correspond- 
ing decay rates are Ai + Rc Ci and Ai — Re Ci, respectively, 
with the dipole-dipole coupling parameters Ci given explicitly 
in reference lEoll . The quasi-steady state population of |e2) is 
unity, so the transition rate is poi = 2^41, in agreement with 
the result of reference L20I . The other transition rates are a bit 
more complicated. For pio one has to take into account the 
decays from |s23) and |a23) to [62), forpi2 the decays from 
|s23) and |a23) to |si3) and from |si2) and \ai2) to \g), and 
forp2i the decays from |e3) to IS23) and |a23) and from |si3) 
and |ai3) to |si2) and |ai2). Multiplying for each decay the 
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lea) 



:^ I «23)-^ ^ [S23 ) 



-|e2) 



The conditional Hamiltonian, without detuning, and the reset 
state in this case are given by 



3 4 



h.c. 



(21) 



and 



FIG. 3: Level configuration of two D systems in the Dicke basis. 
Transitions with rate A2 ±Re C2 (dotted arrows) and transitions with 
rate Ai ± Re Ci (dashed arrows). Fast transitions (with A3 ± Re C3) 
and line shifts due to detuning and to Im d are omitted. 



decay rate by the steady state population of the initial level 
and adding up the different contributions then yields the same 
results for the transition rates as obtained by the Bloch equa- 
tion approach in reference 1 20]. The same is also true for three 
dipole-interacting D systems I22ll . 



IV. THREE DIPOLE-INTERACTING FOUR-LEVEL 
SYSTEMS 

An application of the simplified method to the four-level 
system describing Ba+ is also possible. As depicted in Figure 
[2 the transition from a bright to a dark period is a two step 
process for this system, first an excitation to level |4) by in- 
coherent light with the rate W and then a decay to level |2) 
with the Einstein coefficient A2. Instead of a single Einstein 
coefficient one therefore has to use the product of the incoher- 
ent transition rate W with the branching ratio ^2/(^2 + A4) 
for a decay from state |4) to state |2) for this transition. Then 
everything works as in the case of the D systems and one con- 
firms the results for a single four-level system aheady known 
from the Bloch equation approach |20|. 

Consequently it is also possible to obtain the transition rates 
for three four-level systems which would be rather involved 
to do with the Bloch equation approach. The Bloch equations 
can be written in the compact form 



P 



ffcondP - pHlJ + (P) + ^(P) (19) 



= {^o + ^i{Ai,W)}p, 

where i^wip) describes the incoherent driving as in reference 
1E5I1 and is given by 

3 

{P) = WY, {S^4PSu + S^aPS^,) , (20) 



with 



S+ = |2),,(1|, S+ = |4),,(2|, S+ = |3),.(1| 
5+ = |4),,(l|, and 



np) = T.T.^^s^,p^^j 

i=i j=i 

3 4 



k,l — l -)■ — 1 
k<l 



(22) 



where 



3A 



1 ia"' 



^kl 



(1 - cos^ 6*, 



kl 



(23) 



,0)2 



,0)3 



(l-3cos2 9ki) 



is the coupling parameter which describes the dipole-dipole 
interaction between atom k and atom I for the transition con- 
nected with the Einstein coefficient Aj, with 6ki being the 
angle between the dipole moments and the line connecting 

the atoms. The dimensionless parameter a^:';' = 2T:rki/Xj is 
given by the inter-atomic distance rki multiplied by the wave 
number 27r/Aj of this transition. In order to get a maximal ef- 
fect of the dipole-dipole interaction we assume as in |22] that 
the atoms form an equilateral triangle (i.e. rki ~ r) and that 
()ki — tt/2. Then C^^ -* becomes the Cj of reference | 2fll. 

The quasi-steady states are already known from the calcu- 
lations for three three-level systems. As in reference f23\ one 
can use a symmetrized basis analoguous to the Dicke basis for 
two atoms. This leads to the states 



\Sijk) 



\aijk) — 



1 

+ 
1 

1 



Oijk) 



{m\k) + \j)m) + \kmj) 

mk)\j) + \jmk) + \k)\j)\z)), (24a) 

m)\j)-\jmk)^\k)\m), (24b) 

2\^)\m-\J)\k)\^)-\kmJ) 



2\i)\k)\j)-\m\k)~\k)\m), (24c) 



h,k)^:,{\j)\m)-m)\j) 



\dijk) 



m\k) + \kmz)), 
{2m)\k)-\j)m) 
2m)\j) + \jmk) 



\km\3) 

mm, 



(24d) 



(24e) 
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jk)^^{\j)\m)-m)\j) 
+ \jmk)~\k)\m), 



(24f) 



i < j < k; i,j,k = 1, ... ,4, in the case where all three 
atoms are in different states. For the remaining states one gets 
fori,j = 1, . . . ,4, i 7^ j. 



1 

1 

7E 
1 

71 



m)\3) + \3)\m + \j)\m) 



(25a) 



'mm~\jm^)~\j)\m) (25b) 



{\m\^)-\j)\m) 



(25c) 



if two atoms are in the same state and 

|5) = |1)|1)|1), \e,) = m)\i) for z = 2,3,4 (26) 

if all three atoms are in the same state. The quasi-steady states 
for intensity periods Jq to I2 are, by symmetry, given by 



P.ss,o = |e2)(e2| 

Ass,i-^{Psf® 12)22(21® 12)33(21 

,1D 



(27a) 
(27b) 



|2)n(2| ® pif ® 12)33(21 + |2)ii(2| ® 12)22(21 ® pI^} 



1 ^ 

P.ss,2 = :r5]pi°2®|2)«(2|, 



(27c) 



i=l 



where p'f is the quasi-steady state of one D system in the 
{|1), |3)} subspace and p^fj is the quasi-steady state in the 
subspace corresponding to double intensity of two D systems. 
The state pss.s is rather complicated. Therefore only the pop- 
ulations of the relevant levels will be given, i.e. 



I 27.521 



I 

O 



27.52 



27.519- 




4 , , 6 



FIG. 4: Double jump rate nui and for three dipole-interacting four- 
level systems with experimental parameter values of reference IT^ . 
Dotted line: independent systems. Cooperative effects are less than 
l%o for distances larger than A3. 



Now the procedure is the same as described in the previous 
section for two D systems and one obtains 



Poi = 3^1 P12 = 2Ai 



P23 = 



and 



Pw 



P21 



A2W{Al 



{A2 + AiM 
2A2W 
Ao + Ai 



^3 



(29a) 



(29b) 



P32 = 



3A2W 
A2 + A4, 



2^1 



Al 



2ni 



2ReC3 



4RcC3 



[A 



2nlf 



(29c) 



^3^3 



2nif 



(29d) 



(ff|Pss,3|5) 



1 r 



{{Al + nl)[{Al + nlf + 2.AlB] 



+2A3[|C3nA3 + C3p + s2]} 



(28a) 



(sii3|As.3|sii3) = -^[{Al + nl){iAl + nl) + iAlB] 

(28b) 

(^113 |Pss,3 1^113) = (Cll3|pss,3|cil3) = ^(A^+l^^) (28c) 



(Sl33|Pss.3|Sl33; = ^(3^3 + ^l) (28d) 
(e3|Pss,3|e3) = (^133|Pss,3|^133) = (ci33 |Pss,3 |ci33) = 

(28e) 



with 



and 



TV = {{Al + 2nl)[{Al + 2nlf+iAlB] 
+2A3[|C3nA3 + C3p + i?2]} 



B 



|C3p + 2A3ReC3. 



as transition rates up to first order in C3. The exact results 
including detuning are given in the appendix. The approx- 
imations to first order in C3 have the same structure as for 
three dipole-interacting three-level systems given in reference 
i22l . Basically this means an increase of cooperative effects 
by a factor of two compared to two atoms. In terms of these 
transition rates the double and triple jump rate, i.e. the rate 
of two or three subsequent jumps within a short time window 
Tw, are then given by |20] 



P0lP2lP32(P01 +P12) 



P2lP32(P01 +P10) +P0lPl2(P23 +P32) 



Tw (30) 



and 



TT-TJ = 2 



POIPIOPI2P2IP23P32 



P2lP32(P01 +J5lo) +PmPl2{P23 +P3.2) 



rp2 



(31) 



In Fig. Ela plot of ?idj for the experimental parameter values 
of reference fl9] is shown. The effects of the dipole-dipole 
interaction are negligibly small in particular for experimental 
distances of about ten times the wavelength A3 of the strong 
transition. Without detuning A3, maximal cooperative effects 
are obtained for 1^3 



rV\/5- IA3. This case is shown in 
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FIG. 5: Triple jump rate nxj for three dipole-interacting four-level 
systems. Parameter values as in Fig. |4| except for A3 — and 
^3 = ^\/V5 — IA3 for maximal effects. Dotted line: indepen- 
dent systems, dashed line: up to first order. Cooperative effects are 
less than 5% for distances larger than A3. 

Fig.|5]for the triple jump rate rijj. For inter-atomic distances 
larger than one wavelength A3 of the strong transition cooper- 
ative effects are less than 5% and again rapidly decreasing for 
larger distances. For non-zero detuning the maximally achiev- 
able effects have about the same value. Also one has to bear 
in mind that, as in reference |22], this result has to be seen as 
an upper limit for all possible configurations in the trap. Large 
cooperative effects, i.e. enhancements of the double and triple 
jump rate by several orders of magnitude, can therefore not 
be explained by the dipole-dipole interaction. Furthermore 
one sees that the first order results of equation ( I29> are a very 
good approximation to the exact transition rates given in the 
appendix. 

V. V SYSTEM AND SIMILAR LEVEL SCHEMES 

From the previous results the question arises whether the 
method presented here is also applicable to level systems like 
the V system, i.e. systems in which the transition between 
different light and dark periods results from a coherent ex- 
citation. It turns out that for these systems the situation is 
much more complicated. For a single V system for example, 
{il.2) contains coherences between the ground state 1 1) and 
the metastable state |2). Therefore ^ipi has no component in 
the subspace of eigenstates of for eigenvalue zero and the 
state at time + At in the Bloch equation approach is thus 
given by O 

pito + At, po,») = Pss,, + (e - ^o^^■^lP,s,^■ (32) 

An explicit evaluation of this expression for a single V sys- 
tem starting with pi not only leads to terms proportional to 
the quasi-steady state population of the ground state but also 
to terms proportional to the quasi-steady state coherence be- 
tween ground state and excited state. 



The situation gets even more involved for dipole-interacting 
V systems. Here the term (e — -Sfo)^^ gives rise to additional 
factors which depend in a very complicated way on C3. This 
is in contrast to the D and the four-level system, for which the 
C3 dependence in the transition rates is solely due to the C3 
dependence of the quasi-steady states. The physical reason 
for this is that the efficiency of the laser driving is influenced 
by the dipole interaction, for example via additional detun- 
ings. Therefore the mechanism of jumps in the light intensity 
based on laser driven transitions is much more complex than 
for jumps based on spontaneous decay and incoherent driv- 
ing so that the method outlined above is applicable only in the 
latter case. 



VI. CONCLUSIONS 

In this paper we have presented a simplified approach for 
the calculation of the transition rates between periods of dif- 
ferent intensity of a system of dipole-dipole interacting atoms 
which show macroscopic quantum jumps in their fluores- 
cence. This method works for atoms with level configurations 
in which the transition between the different intensity periods 
is based on incoherent processes. Results previously obtained 
with other methods are recovered by the new approach. 

In addition, the new method has allowed the calculation 
of the transition rates for three interacting four-level systems 
modeling the the relevant level structure of Ba+ ions. This 
allows a direct comparison with the experiment of reference 
lll9ll . This experiment reported an enhancement of the double 
and triple jump rate by several orders of magnitude and this 
was explained through cooperative effects due to the dipole- 
dipole interaction between the ions. With the present results it 
is seen that this cannot be the explanation for the reported en- 
hancement. Cooperative effects can indeed be found for this 
system but they are much smaller, namely only maximal 5% 
of the values for independent atoms. For the parameter values 
of the experiment they are practically absent. 



APPENDIX A: EXACT TRANSITION RATES INCLUDING 
DETUNING 

As mentioned above the transition rates between the differ- 
ent intensity periods can be calculated exactly in C3 and with 
inclusion of a possible detuning of the laser A3 with respect to 
the corresponding atomic transition. The result for the down- 
ward rates is 

A,W{Al + nl + AAl) 
{A2 + A4)[Al + 2nl + iAl] ^ ' 
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P21 



2A^W {A 



-4Ai)(Ai + 2f]i + 4Ai) + (A 



4Ai)(|C3p 



2A3ReC3 -4A3lmC3) 



P32 



^2 +v44 
2A2W 

' A2 + Ai 

2.A2W 

' Ao + Aa 



4A2 



^2 + 2nl + 4A2 



2^1 



4A2)2 



2RCC3 



-4Ai)2 



h (A2 + 4A2)(|C3P + 2^3ReC3 - 4A3lniC3 
A3f^2(A2+4A2) 



2ni 



-4Ai]3 



4 Im C? 



A3»^(Ag+4A^) 
[^i + 2f^i + 4Ai]3 



0(Cf ). 



-3(^1 



-4Ai)B] 



-2{Al 



-4Ai)[|C3nA3 
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2nl 
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2nl 



AAl 
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4AI) 
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-4AI) 



[|C3|2|A3-2iA3 + C3|2 + S2] 



[A 



2ni 
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- 8 Im C\ 



..nliAl 
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[Al 
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with B = IC3 p + 2A3Re C3 — 4A3lm C3. The upward rates independent of C3 and A3, 
of equation (|29aJ are akeady the exact results since they are 
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